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We construct a family of viable scalar-tensor models of dark energy (DE) which possess a phase 
of late-time acceleration preceded by a standard matter era, while at the same time satisfying the 
local gravity constraints (LGC). The coupling Q between the scalar field and the non-relativistic 
matter in the Einstein frame is assumed to be constant in our scenario, which is a generalization of 
f(R) gravity theories corresponding to the coupling Q = — l/y/6. We find that these models can 
be made compatible with local gravity constraints even when \Q\ is of the order of unity through a 
5_l ' chameleon mechanism, if the scalar-field potential is chosen to have a sufficiently large mass in the 

high-curvature regions. We show that these models generally lead to the divergence of the equation 

^5 . of state of DE, which occur at smaller redshifts as the deviation from the ACDM model become more 

significant. We also study the evolution of matter density perturbations and employ them to place 

f ' bounds on the coupling \Q\ as well as model parameters of the field potential from observations of 

the matter power spectrum and the CMB anisotropics. We find that, as long as \Q\ is smaller than 
the order of unity, there exist allowed parameter regions that are consistent with both observational 
and local gravity constraints. 

Q ! PACS numbers: 98.70.Vc 

to: 

I. INTRODUCTION 

T— I ■ 
> ' 

\Q> The origin of dark energy (DE) has persistently posed one of the most serious mysteries in modern cosmology [l|, Q ■ 

The first step toward understanding the nature of DE is to clarify whether it is a simple cosmological constant or 
it originates from other sources that dynamically change in time. The dynamical DE models can be distinguished 
from cosmological constant by studying the variation of the equation of state of DE (= wde) as well as the evolution 
of density perturbations. The scalar field models of DE such as quintessence @ and k-essence Q predict a wide 
variety of variations in u>de, but still the current observational data are not sufficient to rule out such models unless 
the equation of state shows a peculiar evolution. Moreover, the scalar field is required to have a light mass m$ 
^\ ' comparable to the present Hubble parameter (m^ ~ 10~ 33 eV), in order to give rise to an accelerated expansion. 
^ . This requirement is generally difficult to reconcile with fifth-force experiments unless there exists some mechanism by 
which the interaction range of the scalar-field mediated force can be made shorter. 

There exists another class of dynamical DE models that modify Einstein gravity. The simplest models that belong 
to this class are those that are based on the so called f(R) gravity theories in which the Lagrangian density / is 
a function of the Ricci scalar R. ft is well known that theories of the type f(R) = R + aR 2 can give rise to an 
inflationary expansion in the early universe because of the dominance of the aR 2 term pj. In the context of DE, 
the model f(R) = R — /i 2 ( n+1 ) /R" (n > 0) was proposed to explain the late-time accelerated expansion due to the 
dominance of the term (i 2 ( n+1 * > / R n Q (see also Refs. J?})- It was found, however, that this model is plagued by a 
number of problems such as the instability of matter perturbations [8| as well as the absence of a matter-dominated 
epoch @. 

In the past few years there has been a burst of activity in the search for viable f(R) DE models [Tol. [Til H3. [13. 
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In Ref. [ll|] the conditions for the cosmological viabilities of f(R) DE models (having a matter era 
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followed by an accelerated epoch) were derived without specifying the form of f(R). A number of general conditions 
are required on general viability and stability grounds. For the existence of a prolonged saddle matter era, the quantity 
rn = Rf,Rii/f,R needs to be positive and close to 0. To avoid anti-gravity, f t R is required to be positive in regions 
R > Ri, where R\ (> 0) is a Ricci scalar at a de-Sitter attractor responsible for the accelerated expansion. Also, to 
ensure that density perturbations do not exhibit violent instabilities, we require f^RR > o dm. The conditions 
f,RR > and f t R > (for R > Ri) have been shown to also ensure the absence of ghosts and tachyons [l5j |. 

The local gravity constraints (LGC) should also be satisfied for the viability of f(R) models [20]. The f(R) gravity 
in the metric formalism is equivalent to scalar-tensor theory with no scalar kinetic term, namely, the Brans-Dicke 
model with a potential and wbd = [21]. If the mass of the scalar- field degree of freedom always remains as light as 
the present Hubble parameter H , one can not satisfy the LGC due to the appearance of the long-ranged fifth force. 
It is possible to design the field potential so that the mass of the field is heavy in a large-curvature region where local 
gravity experiments are carried out. Then the interaction range of the fifth force becomes short in such a high-density 
region, which allows the possibility of the models being compatible with LGC. 

In fact a number of viable models based on f(R) theories have been proposed [TH, [H, [13, EH that can satisfy both 
the cosmological and local gravity constraints discussed above. In the high-density region (R 3> R c ) these models 
have asymptotic behavior f(R) ~ R — fJ,R c [l — (R/ R c )~ 2n ] (/i > 0, R c > 0, n > 0), where R c is of the order of the 
present Ricci scalar. Inside a spherically symmetric body with an energy density p m , the field acquires a minimum 
at R ~ p m with a mass much heavier than Hq. In this case the body has a thin-shell inside it so that the effective 
coupling between the field and the matter decreases through the so-called chameleon mechanism [12J, |22|, |23|, 12J, l25j . 
The bounds on the model parameters of such models derived from solar-system and equivalence principle constraints 
are given by n > 0.5 and n > 0.9, respectively [25[. For viable f(R) models there are also a number of interesting 
observational signatures such as the divergence of the equation of state of DE [l8|, [26J and the peculiar evolution of 
matter perturbations [HI, [lj| [l8|, [lj| . This is useful to distinguish f(R) gravity models from the ACDM model. 

In the Einstein frame the f(R) gravity corresponds to a constant coupling Q = — 1 /\/6 between dark energy and the 
non-relativistic fluid @ (see Eq. (fTU)l for the definition of Q) . Basically, this is equivalent to the coupled quintessence 
scenario [27j with a specific coupling. Our aim in this paper is to generalize the analysis to scalar-tensor theories 
with the action (|13| in which the coupling Q is an arbitrary constant. We regard the Jordan frame as a physical 
one in which the usual matter conservation law holds. The dark energy dynamics in scalar-tensor theories has been 
investigated in many papers [H, [29|, [3(1 HH, HH after the pioneering works of Refs. [H, |34|. If the mass of the field <fi is 
always of the order of H , the solar-system constraint wbd > 4.0 x 10 4 [HI] gives the bound |Q| < 2.5 x 10 -3 . Previous 
studies dealing with the compatibility of the scalar-tensor DE models with LGC have restricted their analysis to this 
small coupling region [3(3, [SH- We wish to extend the analysis to the case in which the coupling |<5| is larger than the 
above massless bound. In fact one can design the potential V(4>) so that the mass of the field is sufficiently heavy 
in the high-density region to satisfy LGC through the chameleon mechanism. We shall construct such a viable field 
potential inspired by the case of the f(R) gravity and place experimental bounds on model parameters in terms of 
the function of Q. 

We shall also study the variation of the equation of state for DE and the evolution of density perturbations in 
such scalar-tensor theories. Interestingly, we find that the divergent behavior of u>de is also present as in the case of 
f(R) gravity. We also estimate the growth rate of matter perturbations and show that the non-standard evolution of 
perturbations manifests itself from a certain epoch (depending upon model parameters) during the matter era. This 
is useful to place constraints on model parameters using the data of large scale structure and CMB. 

This paper is organized as follows. In Sec. |TT] we consider a class of scalar-tensor theories with constant coupling 
Q. In Sec. IIIII we study the background cosmological dynamics and consider the case of constant as well as varying 
A (the slope of the potential in the physical frame). In this section we also introduce a family of potentials which 
are natural generalizations of a viable family of models in f(R) gravity. In Sec. IIVI we discuss the LGC under the 
chameleon mechanism and place experimental bounds on parameters of viable scalar-tensor models using solar-system 
and equivalence principle constraints. In Sec. [V] we study the evolution of the equation of state of DE and show that 
the divergence of wde previously found in f(R) theories is also present in the class of scalar-tensor models considered 
here which are compatible with LGC. In Sec. I VII we discuss the evolution of density perturbations and place constraints 
on the coupling and model parameters employing the predicted difference in the slopes of the power spectra between 
large scale structure and the CMB. Finally we conclude in Sec. I VIII The stability analysis which is crucial to derive 
the background cosmological scenario is briefly summarized in Appendix IA 11 In Appendix IA 21 we also clarify the 
stability condition of a de-Sitter point that appears in the presence of the coupling Q for varying A. 
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II. SCALAR-TENSOR THEORIES 



We start with a class of scalar-tensor theories, which includes the pure f(R) theories as well as the quintessence 
models as special cases, in the form 



d x\f— ~g 



(1) 



Here, / is a general differentiable function of the scalar field ip and the Ricci scalar R, £ is a differentiable function of 
<p and S m is a matter Lagrangian that depends on the metric g^ v and matter fields \I/ m . We also choose units such 
that k 2 = 8irG = 1, and restore the gravitational constant G when it makes the discussion more transparent. 
The action |T]) can be transformed to the so called Einstein frame under the conformal transformation [3t| : 



_ 2Q 



where 



= - In F, 
2 



df_ 
dR 



In the following we shall consider F to be positive in order to ensure that gravity is attractive. 
We shall be considering theories of the type 

f&,R)=F(<p)R-2V(<p), 

for which the conformal factor Q depends upon ip only. Introducing a new scalar field (f> by 



3 (F# 
F 



C_ 

F 



the action in the Einstein frame becomes [3 

s r I d'v 



where a tilde represents quantities in the Einstein frame and 

V 



dip, 



+ S m (g flu F 1 ,* m ). 



U 



F 2 



(2) 



(3) 



(4) 



(5) 



(6) 



(7) 



In f(R) gravity theories without the field (/?, the conformal factor f2 depends only on R. Introducing a new scalar 
field to be 



^=^lnF, 

the action in the Einstein frame is given by (O with the potential [36 

RF-f 



U 



2F 2 



(8) 



(9) 



Hence, the f(R) gravity can be cast in the form of scalar-tensor theories of the type (JlJ with (j4]), by identifying the 
potential in the Jordan frame to be V — (RF — f)/2. 

In order to describe the strength of the coupling between dark energy and a non-relativistic matter, we introduce 
the following quantity 



Q 



2F 



(10) 



From Eq. ([H]) one has F — e 2 */v^ wn i cn shows that the f(R) gravity corresponds to 



Q = -1/V6. 



(11) 
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In what follows we shall study a class of scalar-tensor theories where Q is treated as an arbitrary constant. This class 
includes a wider family of models, including f(R) gravity induced gravity and quintessence models. Using Eqs. ([5]) 
and (|10p we have the following relations 



Then action (TT|) in the Jordan frame together with (HJ) yields 



s = 



-9 



\fr- 1 -{i-qq 2 )f{Vcp) 2 -V 



(12) 



(13) 



In f(R) gravity the kinetic term of the field 4> vanishes with the potential given by V = (RF — f)/2. Note that in 
the limit, Q — > 0, the action (|13[) reduces to the one for a minimally coupled scalar field (f> with a potential V(c/>). 
It is informative to compare (I13|) with the following action 



S = 



-9 



1 D ^BD /„ ^ 



V 



(14) 



which corresponds to Brans-Dicke theory with a potential V. Setting \ = F = e 2C ^, one easily finds that two 
actions are equivalent if the parameter wbd is related with Q via the relation 



3 + 2w B d = 



1 

2Q 2 



(15) 



Under this condition, the theories given by (fT3|) are equivalent to the Brans-Dicke theory with a potential V. 

In the following sections we shall in turn consider the evolution of the background dynamics in homogeneous 
settings, the local gravity constraints and the matter density perturbations. 



III. HOMOGENEOUS COSMOLOGY 



In what follows we shall discuss cosmological dynamics for the action ([13]) in the flat Friedmann-Lemaitre- Robertson- 
Walker (FLRW) spacetime ds 2 = — d£ 2 +a 2 (i)dx 2 , where t is cosmic time and a(t) is the scale factor. As a source of 
the matter action S m , we consider a non-relativistic fluid with energy density p m and a radiation with energy density 
p r ad- Then the evolution equations in the Jordan frame are given by 

3FH 2 = ~(1 - 6Q 2 )F0 2 +V- 3HF + p m + p rad , (16) 

2FH = -(1 - 6Q 2 )F4> 2 -F + HF- Pm ~ ^ rad , (17) 
Pm + 3H Pm = , (18) 

Prad + 4ffp r ad = , (19) 

where H = a/a and a dot represents a derivative with respect to t. 
Taking the time-derivative of Eq. P^|) and using Eq. (| 1 7[) , we obtain 




(1 - 6Q 2 )F | o - ;]IJ<>- —,, | + U,, -!- (}FR . (20) 
where the Ricci scalar is given by 

R = 6{2H 2 + H) . (21) 

We regard the Jordan frame as a physical one, since the usual matter conservation holds in this frame [see Eq. (fT5|) ] . 

In order to study the cosmological dynamics, it is convenient to introduce the following dimensionless phase space 
variables 




' != 7fe' " 2 = ^VaF' X3 = li\l^- (22) 
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Then the constraint equation (| 16|) yields 



a, 



1-(1-6Q 2 ) 2\f%Qx\ 



3FH 2 

We also define the following quantities 

f^rad = x i , ^de = (1 - 6Q 2 )x 2 + x\ + 2V6Qxi 

Eq. (f^5]) then yields the relation f2 m + £l ra( j + fi^E = 1 ■ 
From Eqs. (JTTj) and (J20]) we obtain 



'3 • 



H 



1-6Q 2 



3 + 3x\ - 3x1 + x 2 3 - 6Q 2 x 2 + 2V6Q:ri] + 3Q(\xj - 4Q) , 



-^5- = 3(Aa;? - VQx x ) + 3Q (5 - 6Q 2 )a; 2 + 2yf§Qx 1 - 3x 



2,2 i 
1 I, — 1 



Using these relations, we obtain the following autonomous equations: 
dxi 

— (2Q - \)xix 2 - x 2 — 



y ^ (Ax 2 - Vtei) + [(5 - QQ 2 )x\ + 2V6gxi - 3x 2 + x\ - 1 



dx 2 
dF 



= y/QQ XlX3 - 2x 3 - x 3 -^ , 
where N = In (a) is the number of e-foldings and A is defined by 



Xi 



H 
IP 



(23) 



(24) 



(25) 
(26) 

(27) 
(28) 
(29) 



A= — 



Ya 

V 



(30) 



The exponential potential V{<j>) — Vbe~ A( ^ gives a constant value of A. Generally, however, A is dependent on 0, where 
the field is a function of xi, x 2 and x 3 through the definition of x 2 and Eq. ([23]) . Hence Eqs. (|2"7| - (f2"9")) are closed. 
The effective equation of state is given by 



2 H 
1 ~ 3lP 

+ 



(3 + 3x\ - 3x 2 2 



QQ 2 x\ + 2V6Qx 1 ) - 2Q{\x 2 2 - 4Q) . 



(31) 



In what follows we shall first discuss the case of constant A and then proceed to consider the varying A case. 



A. Constant A 



If A is a constant, one can derive the fixed points of the system by setting the r.h.s. of Eqs. (|2~7| - (f2"9"|) to be zero. In 
the absence of radiation (x 3 — 0), we obtain the following fixed points: 



• (a) <fi matter-dominated era (<^>MDE (27J) 

\/6Q 



(xi,x 2 ) = 



3(2Q 2 -1) : 



3-2Q 2 
3(1 - 2Q 2 ) 2 ' 



4Q 2 



3(1 -2Q 2 ) 



(32) 



(bl) Kinetic point 1 



(xi,x 2 ) 



V6Q + 1 



0, w. 



off 



3-V6Q 
3(1 + V6Q) 



(33) 
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(b2) Kinetic point 2 



(xi,x 2 ) 



V6Q-1 



(c) Scalar-field dominated point 



(31,2:2) 



76(4Q - A) 
6(4Q 2 -QA- 1) 

• (d) Scaling solution 

(X1,X2 



6 - A 2 + 8QA - 16Q 2 



6(4Q 2 - QX - 1) 



l/2\ 



^2A' 



3 + 2QA - 6Q 2 



2A 2 



1/2N 



• (e) de-Sitter point (present for A = AQ) 

(xx,x 2 ) = (0,1), Q r , 



»cfF 



\/6g 



3(1 - V6Q) 



TO = , w cff = - 



20Q 2 - 9QA - 3 + A 2 
3(4Q 2 -QA- 1) 



3 - 12Q 2 + 7QX 
A 2 



Wcff 



2Q 
A 



0, W e ff = -1. 
(0,0,1). 



(34) 



■(35) 



(36) 



(37) 



Note that, when X3 7^ we have a radiation fixed point {x\, X2, X3) 

One can easily confirm that the de-Sitter point exists for A = AQ, by setting <j) = in Eqs. Q16p . ((TTJ) and (|20[) . 
This de-Sitter solution appears in the presence of the coupling Q. Note that this is the special case of the scalar- field 
dominated point (c). 

Now given a value for A, and using the stability conditions for the above fixed points given in the Appendix A, the 
cosmological dynamics can be specified. We shall briefly discuss the cases Q = and Q ^ in turn. 



1. Q = 

When Q = (i.e., F = 1, which corresponds to a standard minimally coupled scalar field), the eigenvalues fix and 
/i2 of the Jacobian matrix for perturbations about the fixed points reduce to those derived in Ref. (37j with 7 = 1 
(see Ref. [38[ for earlier works) . In this case the matter-dominated era corresponds to either the point (a) or (d) . The 
point (a) is a saddle node because Hi = —3/2 and ^2 = 3/2. The point (d) is stable for A 2 > 3, in which case fl m < 1. 
The late-time accelerated expansion (w c b < —1/3) can be realized by using the point (c), whose condition is given 
by A 2 < 2. Under this condition the point (c) is a stable node. Hence, if A 2 < 2, the saddle matter solution (a) is 
followed by the stable accelerated solution (c) [note that in this case f2 m < for the point (d)] . The scaling solution 
(d) can have a matter era for A 2 ^> 1, but in this case the epoch following the matter era is not of an accelerated 
nature. 

2. Q^O 

We next consider the case of non-zero values of Q. Here we do not consider the special case of A = 4Q. If the point (a) 
is responsible for the matter-dominated epoch, we require the condition Q 2 <C 1. We then have fl m ~ 1 + 10Q 2 /3 > 1 
and Wcs — 4Q 2 /3, for the </>MDE. When Q 2 < 1 the scalar-field dominated point (c) yields an accelerated expansion 
provided that —>/2 + 4Q < A < \/2 + AQ 1 . Under these conditions the <^>MDE point is followed by the late-time 
acceleration. It is worth noting that in the case of f(R) gravity (Q — — l/v6) the 0MDE point corresponds to, 
fl m — 2 and u) e g = 1/3. In this case the universe in the matter era prior to late-time acceleration evolves as a oc i 1 / 2 , 
which is different from the evolution in the standard matter dominated epoch Q . 

We note that the scaling solution (d) can give rise to the equation of state, w e g ~ for \Q\ <§; |A|. In this case, 
however, the condition w e g < —1/3 for the point (c) gives A 2 < 2. Then the energy fraction of the pressureless 
matter for the point (d) does not satisfy the condition £l m ~ 1. In summary the viable cosmological trajectory 
corresponds to the sequence from the 0MDE to the scalar-field dominated point (c) under the conditions, Q 2 < 1 
and -V2 + 4Q < A < V2 + AQ. 



1 Note that under the condition Q 2 <C 1 and in the case where the dynamics is in the accelerated epoch, the condition \QX\ < 1 is also 
satisfied. 



7 



B. Varying A 

When the time-scale of the variation of A is smaller than that of the cosmic expansion, the fixed points derived 
above in the case of constant A can be regarded as the "instantaneous" fixed points [39| ■ We shall briefly consider the 
cases of Q = and Q ^ in turn. 



1. Q = 

We begin with a brief discussion of the Q — case. If the condition A 2 < 2 is satisfied throughout the cosmic 
evolution, the cosmological trajectory is similar to the constant A case discussed above except for the fact that the 
fixed points are regarded as the "instantaneous" ones. In this case the saddle matter solution (a) is followed by the 
accelerated point (c). 

When A 2 ^ 1 the scaling solution (d) is stable with il m ~ 1. Hence the cosmological trajectory during the matter 
era chooses the scaling solution (d) rather than the saddle point (a). If |A| decreases at late-times, such that it satisfies 
the acceleration condition A 2 < 2, the trajectory stops following the solution represented by the matter point (d) to 
follow the scalar-field dominated point (c) 2 . A representative model of this type is provided by the double exponential 
potential, V{4>) = Vb(e~ Al0 + e"^), with, A 2 > 1 and A| < 2 [loj]. The assisted quintessence models in Ref. [H| 
also lead to a similar cosmological evolution. 



2. Q^O 

We shall now proceed to consider the case of non-zero Q. If |A| is initially much larger than unity and decreases 
with time, it happens that the solutions finally approach the de-Sitter solution (e) with A = AQ. As we shall show 
in the Appendix IA 21 the de-Sitter point (e) is in fact stable even for the variable A case, if the potential satisfies the 
condition Q(dX/dF) > or d\/d(f> < at the de-Sitter point. 

In the context of f(R) gravity, it has been shown that the model, 



f(R) = R- f iR c [1 - (R/R c )- 2n ] {fi > 0, R c > 0, n > 0) , 



(38) 



is a good example which can be consistent with cosmological and local gravity constraints 18]. Note that the models 
proposed by Hu & Sawicki [l(| and Starobinsky [l5| reduce to this form of f(R) in the high-curvature region (R^> R c ). 

In this model the field (p is related to the Ricci scalar R via the relation e 2 ^/^ = 1 — 2nn(R/ ' R c )~( 2n+1 \ Hence, the 
potential V — (FR — /)/2 can be expressed in terms of the field 4> as 



v(<j>) = 

The parameter A is then given by 



^R c 



1 



2n + 1 



(2n/i) 2 "/( 2n+1 ) 



1 



A = - 



An 



2n + l 



(2n/i) 2 "/( 2n+1 ) 



-i -2n/(2ra+l) 



(l - e 2 ^) 



-l/(2n+l) 



(39) 



(40) 



In the deep matter-dominated epoch in which the condition R/R c 3> 1 is satisfied, the field <p is very close to zero. 
For n and \x of the order of unity, |A| is much larger than unity during this stage. Hence the matter era is realized by 
the instantaneous fixed point (d). As R/R c gets smaller, |A| decreases to the order of unity. If the solutions reach the 
point A = AQ = — 4/Vt5 and satisfy the stability condition dX/dF < the final attractor corresponds to the de-Sitter 
fixed point (e). 

For the theories with general couplings Q, let us consider the following scalar-field potential 



V{<j>) = V [1- C(l - e- 2 Wf] (V > 0, C> 0, < p < 1) , 



(41) 



2 Note that the de-Sitter solution (e) exists only for A = 0, i.e., for the case of cosmological constant (V =const.). 
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as a natural generalization of Eq. (|39| . The slope of the potential is given by 

_ 2CpQe-^(l-e-^r- 1 

1 _ C(l _ e-2Q0 )p ■ W 

When Q > 0, the potential energy decreases from Vq as <f> increases from 0. On the other hand, if Q < 0, the potential 
energy decreases from Vo as <f> decreases from 0. In both cases we have V((j>) — > Vq(1 — C) in the limits <f> — > oo (for 
Q > 0) and — > -oo (for Q < 0). 

In the model (|4Tj) the field is stuck around the value = during the deep radiation and matter epochs. In these 
epochs one has R ~ p m /F from Eqs. p^|) . (fTT)) and (f2"Tj) by noting that Vb is negligibly small compared to p m or p ra d- 
Using Eq. (|20[) . we obtain the relation + Qp m ~ 0. Hence, in the high-curvature region the field evolves along 
the instantaneous minimum given by 

We stress here that a range of minima appears depending upon the large energy density p m of the non-relativistic 
matter. As long as the condition p m ^> VopC is satisfied, we have |0 m | <C 1 from Eq. (|43|) . 

Since from Eq. (|42|) |A| 3> 1 for field values around = 0, the instantaneous fixed point (d) can represent the 
matter-dominated epoch provided that \Q\ <C |A|. The deviation from Einstein gravity manifests itself when the field 
begins to evolve towards the end of the matter era. The variable F = e _2( ^ decreases in time irrespective of the sign 
of the coupling strength and hence < F < 1. This decrease of F is crucial to the divergent behavior of the equation 
of state of DE, as we will see in Sec. IVl 

The de-Sitter solution corresponds to A = 4Q, i.e., 

C = (44) 

(l-fi)f-i [2+(p-2)fi]' 1 ' 

where F\ is the value of F at the point (e). Provided that the solution of this equation exists in the region < Fi < 1, 
for given values of C and p, the de-Sitter point exists. From Eq. (|42[) we obtain 

dA _ AC P Q 2 F{\ - F)p- 2 [1 - V F - C(l - F)p] 

A4>~ [1 - C(l - F)p} 2 ' ( ' 

When < C < 1, one can easily show that the function g(F) = 1 — pF — C(l — F) p is positive in the region 
< F < 1 giving dX/dcj) < 0. Hence, the conditions for a stable de-Sitter point is automatically satisfied. In this case 
the solutions approach the de-Sitter attractor after the end of the matter era. 

When C > 1, the function g(F) becomes negative for values of F that are smaller than the critical value F c (< 1). 
The de-Sitter point (e) is stable under the condition 1 — pFi > C(l — Fi) p . Using Eq. (|44f we find that this stability 
condition translates to 

Fi>^—- (46) 
2 -P 

If this condition is violated, the solutions choose another stable fixed point as an attractor. In f(R) gravity, for 
example, the solutions can reach the stable accelerated point (d) characterized b, m = — r— 1 and (v3— l)/2 < m < 1 
[ll|, where m = Rf,R R /f lR and r = -Rf iR /f. 

In summary, when < C < 1, the matter point (d) can be followed by the stable de-Sitter solution (e) for the 
model (|4Tj) . In Fig. [T] we plot the evolution of ^de, ^m, ^rad and w e s for Q = 0.01, p = 0.2 and C = 0.7. Beginning 
from the epoch of matter-radiation equality, the solutions first dwell around the matter point (d) with w e s — and 
finally approach the de-Sitter attractor (e) with w e s — — 1. We have also numerically confirmed that A is initially 
much larger than unity and eventually approaches the value A = AQ. 



IV. LOCAL GRAVITY CONSTRAINTS 



In this section we shall study the local gravity constraints (LGC) for the scalar-tensor theories given by the action 
(fl~3|) . In the absence of the potential V(<f>) the Brans-Dicke parameter ljbd is constrained to be wbd > 4.0 x 10 4 
from solar-system experiments [35j. Note that this bound also applies to the case of a nearly massless field with the 
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FIG. 1: The evolution of 57de, fim, f^rad and w e ff for the model (|4ip with parameters Q = 0.01, p = 0.2 and C = 0.7 and 
initial conditions x\ = 0, 2:2 = 2.27 x 10~ 7 , xz = 0.7 and xa — 1 = —5.0 x 10 -13 . 

potential V(</>) in which the Yukawa correction e~ Mr is close to unity (where M is the scalar field mass and r is an 
interaction length). Using the bound cjbd > 4.0 x 10 in Eq. (fT5|) , we find 

\Q\ < 2.5 x 10 -3 (for the massless case). (47) 

This is a strong constraint under which the cosmological evolution for such theories is difficult to be distinguished 
from the Q = case. 

Let us then consider the case in which the mass M of the field <f> is sufficiently heavy so that the interaction range 
of the field (~ 1/M) becomes short so as to satisfy LGC. In the context of f(R) gravity (Q — — 1/\/6) it was in fact 
shown that the LGC can be satisfied by constructing models in which M is large enough in a high-dense region where 
local gravity experiments are carried out [TBI, [l(| EE EH, Eil • In these models the mass M tends to become lighter 
with the decrease of the Ricci scalar R towards the present epoch. In what follows, we shall construct viable models, 
based on scalar-tensor theories whose couplings Q are of order unity, which are consistent with LGC. 

A. Chameleon mechanisms 

In Refs. [H, [IB, HH HH| ^ was explicitly shown that in f(R) gravity a spherically symmetric body forms a thin- 
shell inside the body through a chameleon mechanism [22I [23 ] . Generally this happens in a non-linear regime where 
the mass M of a scalar-field degree of freedom is heavy so that the usual linear analysis based on the inequality, 
\SR\ <C |-R^|, is invalid (where SR is a perturbation about a background value R^). In what follows we shall 
briefly review the chameleon mechanism for the theory given in Eq. (|13[) and then place constraints on viable models 
consistent with LGC. 

Let us consider the Einstein frame action ((6|). The variation of this action with respect to <p leads to the following 
equation of motion 

n<t>-U, (t , = -Qf, (48) 

where, T = e A ^T and T — g^T^, with being energy momentum tensor of the matter in the Jordan frame. We 
take a spherically symmetric spacetime with a radius f from the center of symmetry In this setup Eq. (|48[) becomes 
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where 



t/eff($ = U{<t>) 



Orb * 

p 



(50) 

p, which is conserved 



Here, p* is related with the energy density p = —T in the Jordan frame via the relation p* = e 3 ®^ 
in the Einstein frame [23[ (i.e., p*r 3 =constant). 

We consider a configuration in which the spherically symmetric body has a constant density p* — p* A inside the 
body (f < f c ) and that the energy density outside the body (f > f c ) is given by p* = p* B (< Then the mass of 
this body is given by M c = {Air/2>)r^.p A = (4n/3)r 3 p A . Let us denote the field value at the minimum of the effective 
potential U c s{4>) corresponding to the density p* A (p* B ) by <p A (0s)- That is, they are given by 



V, 



+ Qe^ A p* A =Q, and U^b) + Qe Q * B p* B = 



(51) 



respectively. Under the condition p* A 3> p* B the mass squared m A = U^^a) is much larger than m B = £/" ff (0s). 

In solving for a static spherically symmetric field configuration, we impose the boundary conditions d<f>(r = 0)/df = 
and <p(f — > oo) = <pB, so that the field <p is non-singular at f = and that the force on a test body vanishes at a 
sufficiently large distance. Then <f> starts to roll down the potential where the term dU c s/d(f> becomes important. 

If the field value at the center <j>(r = 0) is close enough to the equilibrium value 4>a with \<fi(r = 0) — 4>a\ -C \4>a\, 
the thin-shell solution is realized [23|]. In this case the field does not move away from <fi(r = 0) practically up to a 
radius f\ which satisfies Ar c /f c = (f c — f\)/f c < 1. At f — fi, the field starts to roll down the potential and we find, 
\U,rj>{4>)\ ^ \Q e ®^P*A\ f° r ?i < r < r c . Under the condition \Q4>a\ <C 1 (as we will confirm later) the r.h.s. of Eq. 
is approximately given by dJ7 c g/d0 ~ Qp\. Outside the body (r > r c ) the gradient energies on the l.h.s. of Eq. 
become important because the energy density drops down from p* A to p* B . Taking into account the mass term ms of 
the effective potential U c s, one has dU e g/d(j) = m B ((f> — 4>b) on the r.h.s. of Eq. (f45|). 

We match the solutions of Eq. (|49p at f — f c with the boundary conditions <f> = 4>a, d<f>/dr = at r = f\ and 
4>{f — > oo) = 4>b- Then the following solution is obtained in the region r > r c (23l. [25|. 



(r) 



where 



QMc 

47T 



-mB(f—f c ) 



(52) 



1 - 



- 4>A 



87rf c 



GM C 



(53) 



In deriving the relation J52[) . we assumed the condition m B f c <C 1. Since we are in the thin-shell regime, we obtain 
the following relation from Eq. (|53p . 

X (54) 



Ar c ^ (pB -4 
Then the solution outside the body (r > f c ) is given by 

W - — i 1" 0s , 

47r r 



where Q ff = 3Q 

r r . 



(55) 



Thus, when the body has a thin-shell, the effective coupling \Q e s\ in the thin-shell becomes much smaller than unity, 
even if \Q\ itself is of the order of 1. 

If the field value at f = is not close to 4>a (i.e., \<j>(f = 0) — 4>a\ ~ I0a|), the field rapidly rolls down the potential 
at fi ~ 0. Setting f\ = in Eq. (|52p . we obtain the solution ([53)) with Q c ff replaced by Q. This is the thick-shell 
regime in which the effective coupling is not small as to satisfy the LGC. 

The presence of the fifth force interaction mediated by the field <f> leads to a modification to the spherically symmetric 
metric. Under the weak field approximation, the spherically symmetric metric in the Jordan frame is given by [l2l.l25T] 



ds 2 = - [ 1 - ] dt 2 



1 



2jG cS M c 



dr 2 +r 2 (d9 



2 • sin 2 



r J \ r 

where the effective gravitational "constant" G c s and the post-Newtonian parameter 7 are given by 



G c ff ~ G 



1 



Q cB e 



7 : 



1 + (\/6Qoff/3)(l + m B r)e- m ^ r - r ^ 
1 - {VGQ c s/3)e- m B(r-r c ) 



(56) 



(57) 
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Note that we have used the approximation f ~ r that is valid in the region \Q<j>\ <C 1. 

Provided that the condition m^r < 1 holds in an environment in which local gravity experiments are carried out, 
we have 7 ~ (1 + v / 6<5eff/3)/(l — \/6QcS /3). Hence, if |Q c ff| is much smaller than unity through the chameleon 
mechanism, it is possible to satisfy the following severest solar system constraint that comes from a time-delay effect 
of the Cassini tracking [42J : 

|7- 1| < 2.3 x 1CT 5 . (58) 

Using the thin-shell parameter, this bound translates into 

Ar c 4.7 x 1(T 6 
r c \Q\ 

If the body does not have a thin-shell for \Q\ of the order of unity, the condition (|55|) is not satisfied. In f(R) gravity, 
for example, we have 7 = 1/2 for Q c s = Q — —l/y/6. 

B. Solar system constraints 

In f(R) gravity, the models can satisfy LGC because the mass M of the field potential is sufficiently 
heavy in the high-density region whose Ricci scalar R is much larger than R c . Since the field mass Ma, inside the 
body is much heavier than that outside the body, most of the volume element within the core does not contribute to 
the field profile at r > r c except for the thin-shell regime around the surface of the body (note that this contribution 
is proportional to e~ Me , where £ is a distance from the volume element to a point outside the body). In the case of 
general couplings Q, the models presented in Eq. (|4ip can be compatible with LGC. Under the condition \Q<j>\ -C 1, 
one has U }( p ~ — 2VoQpC(2Q(j)) p ~ 1 for the potential U = V/F 2 in the Einstein frame. Then from Eq. (|5Tj) we obtain 
the field values at the potential minima inside/outside the body: 



1 {2V pC\~ , 1 (2V pC 



2Q\ pa J 1 2Q V PB 



(60) 



which satisfy \4>a\ <C \4>b\- Note that these are analogous to the field value 4> m derived in Eq. (|43[) in the cosmological 
setting. In order to realize the accelerated expansion at the present epoch, Vq needs to be roughly the same order as 
the square of the present Hubble parameter Hq, so we have Vb ~ Hq ~ po, where po ~ 10~ 29 g/cm 3 is the present 
cosmological density. Note that the baryonic/dark matter density in our galaxy corresponds to ps — 10 -24 g/cm 3 . 
This then shows that the conditions, |Q0a| "C 1 and |Q0_b| <C 1 are in fact satisfied provided that C is not much 
larger than unity. 

The field mass squared M? = d 2 U/dcj) 2 at <fi — <pA is approximately given by 

2-P 



'W' oKyw ^lgJ v °- (61) 

This means that M^c^a) can be much larger than H Q due to the condition pa 3> Vq. Therefore, while the mass Ma, 
is not different from the order of Hq on cosmological scales, it increases in the regions with a higher energy density. 

Let us place constraints on model parameters by using the solar system bound (|59[) . In so doing, we shall consider 
the case where the solutions finally approach the de-Sitter point (e). Since we have Ar c /r c ~ cf>B /(6Q$ C ) with (j>B 
given in Eq. ([60|) . the bound ([59| translates into 

(2V a pC/p B ) 1/(1 - p) < 1.2 x 10- 10 |Q| , (62) 

where we have used the value $ c = 2.12 x 10 -6 for the Sun. At the de-Sitter point (e), one has 2>F\H 2 = Vq[1 — C(l — 
F\) p ] with C given in Eq. (|44|) . Hence, we obtain the following relation 

V = 3Hf + ^- 2 ^. (63) 



Substituting this in Eq. (|62p we find 



Rx 

(>B 



1/(1-P) 

-10 



(1 - F x ) < 1.2 x Hr iU \Q\ , (64) 
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where R\ = l2Hf is the Ricci scalar at the de-Sitter point. Since the term (1 — F\) is smaller than one half from the 
condition (1461) we obtain the inequality {Ri/ p B ) 1/[1 ~ p) < 2.4 x 1CT 10 |Q|. We assume that R\ is of the order of the 
present cosmological density po = 10 -29 g/cm 3 . Taking the baryonic/dark matter density to be pb = 1CP 24 g/cm 3 
outside the Sun we obtain the following bound 

P > 1 - ^ • (65) 

9.6-log 10 \Q\ 

For \Q\ = 10~ 2 and \Q\ = 10 _1 this gives p > 0.57 and p > 0.53 respectively. The above bound corresponds to 
p > 0.50 for the case of f(R) gravity, which translates into the condition n > 0.5 in Eq. ([39]) . This agrees with the 
result found in Ref. 12511. 



C. Equivalence principle constraints 



Let us proceed to the constraints from a possible violation of the equivalence principle (EP). Under the condition 
that the neighborhood of the Earth has a thin-shell, the tightest bound comes from solar system tests of the EP that 
makes use of the free-fall accelerations of the Moon (om ooi1 ) and the Earth (a©) toward the Sun [23|, (25|. The bound 
on the differences between the two accelerations is 42] 



ae\ <10 -i 3i 



(66) 



Since the acceleration induced by a fifth force with the field profile (f>{r) and the effective coupling is given by 
flfifth = \Qcs4>{r)\ we obtain [23| 



GMp 



a© = 



Ar 



GMr. 



A? 



r© 



<I> 2 



0^Moon 



(67) 



where $ Q ~ 2.1 x 10~ 6 , <&© ~ 7.0 x 10 -10 and $Moon — 3.1 x 10~ n , are the gravitational potentials of the Sun, the 
Earth and the Moon, respectively. Note that Ar^/r© is the thin-shell parameter of the Earth. From the bound (|66p . 
this is constrained to be 



< 



3.8 x 10 



-7 



(68) 



Note also that the thin-shell condition for the neighborhood outside the Earth provides the same order of the upper 
bound for, Ar e / r e [25] . 

Taking a similar procedure as in the case of the solar system constraints discussed above (using the value R\ = 10~ 29 
g/cm 3 and pb = 10~ 24 g/cm 3 ), we obtain the following bound: 



p > 1 



13.8 - log 10 |Q| 



(69) 



This is tighter than the bound l[65|) . When \Q\ = 10 2 and \Q\ — 10 1 we have p > 0.68 and p > 0.66, respectively. 
In the case of f(R) gravity the above bound corresponds to p > 0.65 which translates to n > 0.9 for the potential 

(22). 

In summary, the LGC can be satisfied under the condition (|69|) for the potential (|4"Tj) . 



D. General properties for models consistent with LGC 



In this subsection we shall consider the general properties of scalar-tensor theories consistent with LGC, without 
specifying the form of the field potential. In order to satisfy the LGC we require that \<j>B — 4>a\ is much smaller 
than \Q$ C \ from Eq. Since there is a gap between the energy densities inside and outside of the spherically 

symmetric body we have, \<j> B — 4>a\ — \4>b\i which implies \4>b\ "C |Q$ c |- The gravitational potential <I> C is very 
much smaller than unity in settings where local gravity experiments are carried out, hence this yields the constraint 
I 4>b I *C 1. Cosmologically this means that \<p\ is much smaller than unity during matter/radiation epochs. When 
|Q| 3> 1 the condition \4>b\ "C 1 is not necessarily ensured, but those cases are excluded by the constraints from 
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FIG. 2: This illustration describes a field potential V((j>) that is consistent with LGC. For a coupling Q that is positive (negative) 
the potential evolves in the region (f> > (<j) < 0). In the figure (f> m represents the field value during the radiation/matter eras, 
which instantaneously changes in time. The field value <f>^s corresponds to the one at the de-Sitter point. Note that both <f> m 
and 0ds are sustained by the presence of the coupling Q having potential minima characterized by the condition (|70|l . In the 
early stages of the cosmological evolution, the mass M of the field <j> is heavy for consistency with LGC. This mass gradually 
gets smaller as the system approaches the de-Sitter point. 

density perturbations unless the model is very close to the ACDM model (as we shall see later). In the following we 
shall consider the theories with \Q\ < 1. 

In the region \cf)\ <§; 1 (i.e., F ~ 1), the derivative terms are negligible in Eq. (PU|) and the field stays at instantaneous 
minima given by 



in the late radiation-dominated and matter-dominated eras. The condition (|T0[) translates into X/Q — p m /V which 
means that X/Q 3> 1 in the radiation and matter epochs. This is in fact consistent with the condition |w c ff| — 
|2Q/A| <C 1 for the existence of a viable matter point (d). If the de-Sitter point (e) is stable, the solutions finally 
approach the minimum given by (|70[) . i.e., X/Q = 4. 

The sign of A needs to be the same as that of Q in order to realize the above cosmological trajectory. When Q > 0, 
we require A = —V^/V > 0, i.e., V t< j, < 0, which means that the field evolves along the potential toward larger 
positive values from <fr ~ 0. When Q < the field evolves towards smaller negative values from tj> ~ 0. 

We illustrate such potentials in Fig. O Since the ratio X/Q decreases from the radiation/matter epochs to the 
de-Sitter epoch, the derivative dX/d<fi is negative irrespective of the sign of Q. We recall that in this case the stability 
of the de-Sitter point (e) is also ensured. Since dX/d<p = X 2 — V^/V, the mass squared 



is required to be positive to satisfy the condition dX/dcj) < 0. Moreover, the mass M needs to be heavy enough 
in order to satisfy the condition M 2 > X 2 V in radiation/matter epochs. The model (|4"l"T) provides a representative 
example which satisfies all the requirements discussed above. 

It is worth mentioning that for the models that satisfy LGC, the quantity F — e~ 2 Q$ in the matter/radiation eras 
is larger than its value at the de-Sitter point. It is this property which leads to an interesting observational signature 
for the DE equation of state, as we shall see in the next section. 



V i4> + QFR = , 



(70) 



M 2 = , 



(71) 
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V. THE EQUATION OF STATE OF DARK ENERGY 

In scalar-tensor DE models, a meaningful definition of energy density and pressure of DE requires some care. In 
this section, following Ref. [28j . we shall discuss the evolution of the equation of state of DE, which could provide 
comparisons with observations. In the absence of radiation, Eqs. (|16[) and (fTT)) can be written as 

3F Q H 2 = PBE+Pm, (72) 
-2F a H = PDE+PDE + Pm, (73) 

where the subscript "0" represents present values and 

Pde = ^(1-6Q 2 )F4> 2 + V-3HF-3(F-F )H 2 , (74) 
Pde = ^(l-6Q 2 )F<j) 2 -V + F + 2HF+(F-Fo)(3H 2 + 2H), (75) 
which satisfy the usual conservation equation 

Pde + 3i? (pde + Pde) = . (76) 
We define the equation of state of DE to be 

PDE W cS 

wde = = ^ n , (77) 

Pde I - {F/F )il m 

where f2 m and w e g are defined in Eqs. (|23|) and (f3"T|) . respectively. Integrating Eq. (fl"8l) . we obtain 

Pm = 3F Q n^H 2 (l + zf, (78) 

where flm^ is the present energy fraction of the non-relativistic matter and z = ao/a — 1 is the redshift. On using 
Eqs. dZU) and (73} , we find 

3r-(l + z)(dr/dz) 

WDE = , (79) 

3r-3^ } (l + z) 3 

where r — H 2 (z)/Hq. Note that this is the same equation as the one used in Einstein gravity Q. By defining the 
energy density p^ and the pressure pde as given in Eqs. ([74")) and ([75)1 . the resulting DE equation of state wde agrees 
with the usual expression which can be used to confront the models with SNIa observations. 

From Eq. (I77[) we find that u>de becomes singular at the point Q m = Fq/F. This happens for models in which 
F increases from its present value F as we go back in time. From Eq. (|12[) it is clear that F decreases in time for 
Q4> > 0. We note that even when the system crosses the point f2 m = Fq/F physical quantities such as the Hubble 
parameter remain to be continuous. 

The models (|4"Tj) satisfy this condition regardless of the sign of Q, which means that the divergent behavior of u>de 
indeed occurs. We recall that in the context of f{R) gravity (Q = -1/V&) the models f(R) = R - p 2 <- n+1 ')/R n (n > 0) 
correspond to a scalar field potential that decreases toward larger cj), i.e., o > 0|. Hence, the divergence of u>de 
does not occur in such models because of the decrease of F toward the past. 

For the models that satisfy |A| ^> 1 initially such that |A| decreases with time, the solutions are in the regime around 
the instantaneous fixed point (d) during the matter era and finally approach either the scalar-field dominated point 
(c) or the de-Sitter point (e). In Fig. [3] we plot the evolution of u>de for the case Q = 0.1 and C — 0.95, with three 
different values of p. In these cases the final attractor corresponds to the de-Sitter point (e) satisfying the relation 
A = 0.4. During the deep matter era the solutions evolve along the "instantaneous" fixed point (d) with fl m close 
to 1 (because A > 1). After A decreases to the order of unity, the solutions approach the de-Sitter solution (e) with 
Q m = and wue = w cS = -1. 

Figure [3J clearly shows that wde exhibits a divergence at a redshift z c that depends on the values of p. When 
p = 0.3, for example, the divergence occurs around the redshift, z c — 3. For compatibility with LGC we require 
p > 0.53, from solar system constraints, and p > 0.66 from EP constraints, as we showed in the previous section. In 
those cases the critical redshift z c gets larger, which is out of the observational range of current SNIa observations. 
Nevertheless, the equation of state wde shows a peculiar evolution that changes from u>de < — 1 to wde > — 1 at a 
redshift around z c = 0(1). This cosmological boundary crossing, similar to the divergence of wde, is attributed to 



15 




FIG. 3: Figure depicting the evolution of ibde for Q = 0.1 and C = 0.95 with three different values of p (0.3,0.55,0.7). The 
redshift z c at which the divergence of iode occurs decreases for smaller p. 



the fact that F increases as we go back to the past. It is worth noting that this is a common feature among viable 
models that are consistent with LGC, as we have illustrated in the previous section. 

Note that in the limit Q — ► the potential V(4>) approaches a constant value V{4>) — > Vq(1 — C). Hence, the models 
are hardly distinguishable from the ACDM model. In these cases the critical redshift z c also goes to infinity. Thus, 
the effect of modified gravity is more apparent for larger \Q\ and smaller p. In f(R) gravity, for example, the model 
given by Eq. can give rise to the redshift z c as close as a few [18| while satisfying the LGC (p > 0.65). These 
cases are particularly interesting to place tight bounds on model parameters from future high-precision observations. 



VI. MATTER DENSITY PERTURBATIONS 



In this section we discuss the evolution of matter density perturbations and the resulting spectra for scalar-tensor 
theories. Considering as our background spacetime the flat FLRW metric, the perturbed metric including the scalar 
metric perturbations $ and \& in the longitudinal gauge is given by [43j | 

ds 2 = -(1 + 2$)df 2 + a 2 {t)(l - 2V)dx i dx j . (80) 

In the following we shall neglect radiation and consider only the pressureless matter. The components of the energy 
momentum tensor of the pressureless matter are given by 

Tq = -{p m + 5p m ) , Tf = -p m v m ,i , (81) 

where v m is related to the velocity potential V through v m = — (V — b). In the Fourier space matter perturbations 
satisfy the following equations of motion [l^, [53, |44j : 

$ = (82) 
fc 2 

(S Pm / Pm y =3^ - —v, (83) 

where v = av m is a covariant velocity perturbation and k is a comoving wavenumber. We shall introduce the following 
gauge- invariant density contrast: 

S m = — + ZHv . (84) 
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The evolution equation for 8 m is then given by 



2H8 m + ^ = 3B + 6HB . 



(85) 



where B = Hv + ^. 

In Fourier space the scalar metric perturbations in scalar-tensor theories satisfy the following equations \44 



1 

2F 



u6S6+ -(uj^6 2 - F i4> R + 2V c/l )66-3H8F + (SH + 'SH 2 --^] 8F 



(3HF - w0 2 )$ + 3F{H<S> + + 8p„ 



$ + $ = ( + ^i? _ _ F <p + PmV ) ) 



ZH + -^6 ) 86 

LU 



2 



w J ,<j> 2 



2V^ - F i? 
2w 



to ) 2ui 



where u> = (1 — 6Q 2 )F and 



(86) 

(87) 
(88) 



(89) 



SR = 2 



-3(H<S> + *)' - 12H(H$ + *) 



-3i? $-2— * 



(90) 



As long as the mass M defined in Eq. (fTTj) is sufficiently heavy to satisfy the conditions M 2 ^> i? and M 2 > A 2 V 
(in order to ensure dA/d0 < 0), one can approximate ((2V^ — F :< f ) R)/2uj)^ ~ M 2 /uj in Eq. (|89|) . While this quantity 
becomes negative for Q 2 > 1/6 this does not imply that the perturbation 86 exhibits a negative instability. In fact 
we shall illustrate below, that due to the perturbation 8R on the r.h.s. of Eq. (fB"9"]) . the effective mass produced is 
positive. 

Generally, the solution of Eq. (|89f consists of the sum of the matter-induced mode 50i n d sourced by the matter 
perturbation and the oscillating mode S6 osc (scalarons |5|), i.e., 86 = 86^ + S6 osc . The oscillating mode corresponds 
to the solution of Eq. (|8"9")l without the matter perturbation. 

In order to derive the approximate perturbation equations on sub-horizon scales, we use the approximation according 
to which the terms containing k 2 /a 2 , 8p m , 8R and M 2 dominate in Eqs. (f8"6"]) -(|89 p . This method was used in 
Refs. [2l. 1341 14{| in the nearly massless case (M 2 < H 2 ). In the context of f(R) gravity this approximation was shown 
in Ref. [lj| to be extremely accurate even in the massive case (M 2 3> H 2 ) as long as the oscillating degrees of freedom 
do not dominate over the matter-induced mode. 

In order to extract the peculiar features of the matter perturbations in scalar-tensor tensor theories, let us first 
concentrate on the matter induced mode. Under the above-mentioned approximation, we have <5i?i n d — — 2(£; 2 /a 2 )[ , I'-|- 
(I< , ) 0/F)<^i n d] from Eqs. l[8"8"|) and ([M| . where the subscript "ind" represents a matter induced mode. Then from 
Eq. {89]) we find 



2QF 



v ind 



(k 2 /a 2 ){l-2Q 2 )F + M 2 



Using Eq. (SB) and (JHHJ) we obtain 



fc 2 

— * 

or 



8 Pm (k 2 /a 2 )(l-2Q 2 )F + M 2 



2F 



{k 2 /a 2 )F + M 2 



— $ 

o- 



8p m (k 2 /a 2 ){\ + 2Q 2 )F + M 2 



2F 



(91) 



(92) 



(k 2 /a 2 )F + M 2 

which recovers the standard Poisson 



In the limit M 2 /F > k 2 /a 2 one has (k 2 /a 2 )<$> ~ -8p m /2F ~ -4nG8p m , 
equation. In the limit M 2 / F <C k 2 / a 2 one has (fc 2 /a 2 )$ ~ —(8p m /2F)(l + 2Q 2 ), where the effect of the coupling Q 
becomes important. 

From Eq. IjSTjl we find that v is of the order of FH§/p m . Using the fact that (fc 2 /a 2 )$ is of the order of — (l/F)Sp m 
we can estimate that \3Hv/(8p m /p m )\ ~ (aH) 2 /k 2 <C 1. Hence we have 8 m ~ 8p m /p m in Eq. ([5^)1 . Similarly the 
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terms on the r.h.s. of Eq. (|85|) can be neglected relative to those on the l.h.s., which leads to the following equation 
for matter perturbations: 

S m + 2H8 m - 4irG cSPm 5 m ~ , (93) 

where the effective "cosmological" gravitational constant is given by 

_ 1 (fc 2 /« 2 )(l + 2Q 2 )F + M 2 
° ff 8ttF (k 2 /a 2 )F + M 2 ' 1 ' 

We can rewrite Eq. (f9"3")l by using the derivative with respect to N: 

d 2 6 m , (1 3 \ dS m 3 (k 2 /a 2 )(l + 2Q 2 )F + M\ 



dTV 2 V 2 2 / dTV 2 m (k 2 /a 2 )F + M 2 

We also define the effective gravitational potential $ c ff = ($ + ^)/2, which is directly linked with the Integrated- 
Sachs- Wolfe (ISW) effect in CMB and the weak lensing in distant galaxies [l3[ . From Eq. (j9"2")l we obtain the relation 

In order to confront models with weak lensing observations, it is convenient to introduce the anisotropic parameter 
77 defined by, 77 = ($ - [H,|4(|. From Eq. (J92J) we obtain 

„ ±Q\k 2 /a 2 )F 

/_ (fc 2 /a 2 )(l-2Q 2 )F + M 2 ' 1 ' 

which vanishes in the limit M 2 /F 3> k 2 /a 2 but approaches a value 77 — > 4Q 2 / (1 — 2Q 2 ), in the limit M 2 /F <C k 2 /a 2 . 
We also introduce another parameter £ = g(l + rj/2), where q is defined to be (k 2 /a 2 )^ = —(l/2)qp m 8 m . We then 
have S ~ which shows that the effective potential can be written as $ e ff ~ — (a 2 /2k 2 )p m 6 m Y>. Hence, unlike 

the case of the Einstein gravity the weak lensing potential in these scalar-tensor models of gravity are affected by the 
changes of S as well as 5 m . 

During the matter era the field cj> sits at the instantaneous minima characterized by the condition (|70p . This is 

analogous to the situation considered in subsection IV B where for models (|4ip the field value at the potential minimum 

1 

and the mass squared M? are given by Eqs. ([50]) and (fBTj) respectively. Hence, we have the relations <fi oc pm 1 and 

2-P 

M 2 ex M? ex pm p during the matter-dominated epoch. The field 4> can initially be heavy to satisfy the condition 
M 2 /F ^> k 2 /a 2 for the modes relevant to the galaxy power spectrum (0.01ft- Mpc -1 < k < 0.2/iMpc _1 ). Depending 
upon the model parameters and the mode k, the mass squared M 2 can be smaller than k 2 /a 2 during the matter era. 

Next, let us consider the behavior of the oscillating mode. Using Eqs. (|86|) and (188)) under the condition k 2 /a 2 3> H 2 
the gravitational potentials for 5p m = are expressed by 5<fi OBC . Then from Eq. (|90p the perturbation SR corresponding 
to the oscillating mode is given by 

SR osc ~ 6Q \S(t> osc + 3H6<]) osc + ^5<j) osc ) ■ (98) 



Substituting this relation in Eq. (|89|) . we find 

Scj) osc + 3H6cj) osc + + 5ct> OBC ~ , (99) 

which is valid in the regimes M 2 ^> {i?, A 2 U}. Equation (|9"9"|) clearly shows that the effective mass for the oscillating 
mode is positive even for Q 2 > 1/6. 

In the following wc shall confirm that as long as the oscillating mode does not initially dominate over the matter- 
induced mode, it remains subdominant throughout the cosmic history. We shall discuss the two cases: (i) M 2 /F 3> 
k 2 /a 2 and (ii) M 2 /F < k 2 /a 2 , separately. 
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A. The case M 2 /F > k 2 /a 2 

In this regime the matter perturbation equation (|95[) reduces to the standard one in Einstein gravity. During the 
matter era with, w e s ~ and f2 m ~ 1, we have the following solutions: 

S m oc a oc t 2 / 3 , $ e ff = constant . (100) 

2(4-p) 

For the model (|4ip the matter-induced mode of the field perturbation evolves as Sfand oc Sp m /M 2 oc . When 

the frequency uj^ = \Jk 2 /a 2 + M 2 /F changes adiabatically (i.e. |ci^/£J?| <C 1), the WKB solution to Eq. (|99|) is given 
by 




(101) 



For the model (HJ), in the regime M 2 /F 3> k 2 /a 2 , this oscillating mode evolves as S(f> osc oc t 2{1 -^ cos (ct i-pj , where 
c is a constant. 

2 

Now since the background field 4> during the matter era evolves as oc t 1 -p , we find 

5(f)l4> ~ cit 2/3 + c 2 ^ 3 <^ cos (ci - ^) . (102) 

This indicates that the matter-induced mode dominates over the oscillating mode with time. While the solution of 
the oscillating mode in Eq. (|102[) is valid only in the WKB regime (\u! t p/uj 2 1 \ <§; 1), we have checked that 5<p approaches 
a constant value with oscillations at the later stage in which the WKB approximation is violated. Hence, as long as 
the oscillating mode is not overproduced in the early universe, it remains sub-dominant relative to the matter-induced 
mode. Note that this property also holds during the radiation-dominated epoch. 



B. The case M 2 /F < k 2 /a 2 

In this regime the effective gravitational constant (|M|) is given by G e ff = (1 + 2Q 2 )/8irF, which shows that the 
effect of modified gravity becomes important. From Eqs. (|95|) and (|96|) we obtain 

\/25 + 48Q 2 -1 \/25 + 48Q 2 -5 

S m oc t e , $ Gff oc t o , (103) 

which grow faster than the solutions given in Eq. ()100|) . This leads to changes in the matter power spectrum of the 
LSS as well as in the ISW effect in the CMB. 

The field perturbation 5<f> is the sum of the matter-induced mode given in Eq. (|91[) and the oscillating mode 5(j) osc 
given in Eq. (|99|) . Using the WKB solution (jlOip for the latter mode, we have 



V25 + 48Q 2 -5 /o i/o 

54> = Cl t~ § + C2 r 2 l z cos(ct 1/3 ) . (104) 

Since the frequency has a dependence |w^/a;?| ~ H oc l/t, the WKB approximation tends to be accurate at late 
times. Equation (|104|) shows that the matter-induced mode dominates over the oscillating mode with time. 



C. The matter power spectra 



The models (|41j) have a heavy mass M which is much larger than H in the deep matter-dominated epoch, but which 
gradually decreases to become of the order of H around the present epoch. Depending on the modes k, the system 
crosses the point M 2 /F = k 2 /a 2 at t = tk during the matter era. In the context of f(R) gravity this indeed happens 

for the modes relevant to the galaxy power spectrum [TUQjl. Since for the model (T4"Tj) M evolves as M oc < _T ^p 

3(1 -p) 

during the matter era, the time tk has a scale-dependence given by tk oc k i -p . When t < tk, the evolution of S m 
is given by Eq. (|100p . but for t > tk its evolution changes to the form given by (| 1 03[) . 
We define the growth rate of the matter perturbation to be 



_ s m 

3 ~ HS r , 



(105) 
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FIG. 4: The evolution of the growth rate s of matter perturbations in terms of the redshift z for Q = 1.08 and k — 600ao#o 
with three different values of p. For smaller p the critical redshift zt gets larger. The growth rate s reaches a maximum value 
and begins to decrease after the system enters the accelerated epoch. For smaller p the maximum value of s tends to approach 
the analytic value given in Eq. (|106p . 



which is s = 1 in the regime M 2 /F ^> k 2 /a 2 . After the system enters the regime M 2 /F <C k 2 /a 2 during the 
matter-dominated epoch, we have 



g=v ^5+W-l (106) 



During the matter era the mass squared is approximately given by 



r2 , 1 ~P ^fPm\ 1 - p 



M ^W^W^ Q \W v - (107) 

Using the relation p m — 3F S}m^ Hq (1 + z) 3 , we find that the critical redshift Zk at time tk can be estimated as 



Zk 



k 1\ 2(1 p) 2'PpC 1 V 



a H Qj (I-pY-p (3F fi^) 2 -P H o 



4 — p 

- 1 , (108) 



where ao is the present scale factor. The critical redshift increases for larger k/(aoHo). The matter power spectrum, in 
the linear regime, has been observed for the scales 0.01/iMpc -1 < k < 0.2h Mpc" 1 , which corresponds to 30ao-ffo ^ 
k < 600a -ffo- I n Fig- IH we plot the evolution of the growth rate s for the mode k = 600a H n and the coupling 
Q = 1.08 with three different values of p. We find that, in these cases, the critical redshift exists in the region Zk ^ 1 
and that Zk increases for smaller p. When p — 0.7 we have, Zk = 3.9, from Eq. (|108[) . which is consistent with the 
numerical result in Fig. |4j The growth rate s reaches to a maximum value s max and then begins to decrease around 
the end of the matter era. 

McDonald et al. [47| derived the constraint s — 1.46 ±0.49 around the redshift, z = 3, from the measurement of the 
matter power spectrum from the Lyman-a forests. The more recent data reported by Viel and Haehnelt [48| in the 
redshift range 2 < z < 4 show that even the value s = 2 can be allowed in some of the observations. The likelihood 
analysis using these data for the coupled quintessence scenario gives the constraint s < 1.5 [49]. If we use the criterion 
s < 2 for the analytic estimation (|106[) . we obtain the bound Q < 1.08. Figure H] shows that s max is smaller than the 
analytic value s = 2 (which corresponds to Q — 1.08). When p — 0.7, for example, we have that s max = 1.74. For 
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the values of p that are very close to 1, s max can be smaller than 1.5. However these cases are hardly distinguishable 
from the ACDM model. In any case the current observational data on the growth rate s is not enough to place tight 
bounds on Q and p. 

The growth of matter perturbations continues to the time t\ characterized by the condition a = 0. At time t\ the 
matter power spectrum Pg m = (fc 3 /27r 2 )|(5 TO | 2 shows a difference compared to the ACDM model given by 



pACDM y tk 



V25 + 48Q 2 -! 

Ps m {tA) AaVV 6 3 J (l-p)(V25+48Q 2 -5) 



oc k —v . (109) 



The CMB power spectrum is also affected by the non-standard evolution of $ e ff given in Eq. (|103p . This mainly 
happens for low multipoles because of the ISW effect. Since the smaller scale modes in CMB relevant to the galaxy 
power spectrum are hardly affected by this modification, there is a difference between the spectral indices of the 
matter power spectrum and of the CMB spectrum on the scales, k > 0.01/iMpc _1 : 



(i-pK^sTW-s) 

4 — p 

This reproduces the result in the f(R) gravity derived in Ref. [HI]. In Ref. [HI it was further shown that this analytic 
estimation agrees well with numerical results except for large values of p close to unity. This reflects the fact that for 
larger p the redshift z — Zk at time t = tk gets smaller (being of the order of Zk = 0{1\) so the approximations used 
in deriving the solution p03|) . based on w e g = and fi m = 1, break down. In Ref. [18J it was also found that the 
difference An(t ) integrated to the present epoch does not show significant difference compared to (IllOp . 

At present we do not have any observationally significant evidence for the presence of a difference between the 
spectral indices of the CMB and the matter power spectra [5(j. In Fig- E] we plot the constraints coming from the 
criterion, An(t\) < 0.05. If \Q\ is smaller than 0.1, this condition is trivially satisfied. For larger \Q\ the constraints 
on the values of p tend to be stronger. In the f(R) gravity we obtain the bound p > 0.78, which is stronger than 
the constraint coming from the violation of the equivalence principle. If we adopt the criterion An(t\) < 0.03, the 
bound on p becomes tighter: p > 0.87. Meanwhile, if \Q\ is smaller than the order of 0.1, the EP constraint gives the 
tightest bound. If we use the criterion s < 2 for the analytic estimation (|106[) then the coupling \Q\ is bounded from 
above (Q < 1.08). 

In Fig. [5] we show the allowed parameter space consistent with current observational and experimental constraints. 
The constraints coming from the ISW effect in the CMB due to the change in evolution of the gravitational potential 
do not provide tighter bounds compared to those shown in Fig. [SJ 



VII. CONCLUSIONS 



We have considered a class of dark energy models based on scalar-tensor theories given by the action (fl3|) . In 
these theories, expressed in the Einstein frame, the scalar field 4> is coupled to the non-relativistic matter with a 
constant coupling Q. The action (fl~3"]) is equivalent to the Brans-Dicke theory with a field potential V, where the 
Brans-Dicke parameter wbd is related to the coupling Q via the relation 3 + 2ojbd = 1/2Q 2 . These theories include 
the f(R) gravity theories and the quintessence models as special cases where the coupling is given by Q = — l/v6 
(i.e., wbd = 0) and Q = (i.e., wbd — * oo), respectively. 

We began by studying the background cosmological dynamics in a homogeneous and isotropic setting, without 
specifying the field potential V(4>) but under the assumption that the slope of the potential, A = — V^/V, is constant. 
The varying A case can also be studied by treating the fixed points as instantaneous ones. We found that for a range 
of values of the coupling constant \Q\ not much smaller than unity the matter era can be realized by the solution 
corresponding to the point (d) in Eq. ([3l>]) subject to the condition X/Q ^> 1. Interestingly the presence of a non-zero 
coupling Q leads to a de-Sitter solution characterized by the condition + QFR — (i.e., A = 4Q), which can lead 
to late-time acceleration. (The condition for the stability of this de-Sitter solution is given by d\/d(f> < at the fixed 
point.) 

In the absence of the scalar-field potential, solar-system tests constrain the coupling Q to have values in the range 
|Q| < 2.5 x 10~ 3 . The presence of the potential, on the other hand, allows the LGC to be satisfied for larger values 
of |Q|, if the field is sufficiently heavy in the high-curvature region where gravity experiments are carried out. We 
found that even when \Q\ is of the order of 1, a thin-shell can form inside a spherically symmetric body such that the 
effective coupling |Q G ff| defined in Eq. (|55[) becomes much smaller than 1. 

We then considered a family of models given by the scalar-field potentials (|4ip which generalize the corresponding 
potential in the f(R) theory, while at the same time satisfying the LGC for appropriate choices of the parameters. 
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FIG. 5: The allowed region of the parameter space in the (p, Q) plane. We show the bounds coming from the conditions 
An(£/v) < 0.05 and s < 2 as well as the solar-system constraint (|65fl and the EP constraint (|69[) . 



In particular we found that as p approaches unity, the mass of the field (j) becomes larger, thus allowing the LGC to 
be satisfied more easily [see Eq. (fBTj) ]. Using the constraints coming from solar system tests as well as compatibility 
with the equivalence principle, we obtained the bounds p > 1 — 5/(9.6 — log 10 |Q|) and p > 1 — 5/(13.8 — log 10 |Q|), 
respectively. In the f(R) gravity, for example, these constraints correspond to p > 0.50 and p > 0.65 respectively. 

During radiation/matter eras the field <fi needs to be very close to for the compatibility with LGC, which results in, 
F = e~ 2< 3^ ~ 1. Figure [5] summaries the regions of the parameter space in the (p, Q) plane where the corresponding 
potentials lead to models compatible with the LGC. 

For these models we found that the quantity F tends to increase from its present value as we go into the past, which 
results in the equation of state wde of dark energy becoming singular when Q m = Fq/F. This behavior is similar to 
that found for f(R) theories. 

We also studied the evolution of density perturbations for these models in order to place constraints on the coupling 
Q as well as on the parameters of the field potential. In the deep matter era the mass M of the scalar field is sufficiently 
heavy to make these models compatible with LGC, but it gradually gets smaller as the Universe enters the accelerated 
epoch. For the models compatible with the galaxy power spectrum, there exists a "General Relativistic" phase during 
the matter era characterized by the condition M 2 /F ^> k 2 /a 2 . At this stage the matter perturbation 5 m and the 
effective gravitational potential <I> c ff evolve as S m cx t 2 ^ 3 and $ e g = constant, respectively, as in the case of Einstein 
gravity. Around the end of the matter-dominated epoch, the deviation from the Einstein gravity can be seen once 
M 2 /F becomes smaller than k 2 /a 2 . The evolution of perturbations during this "scalar-tensor" regime is given by 
Eqs. (|103[) . Under the criterion s — S m /H5 m < 2 of the growth rate of matter perturbations with the use of the 
analytic estimation (|106[) . we obtain the bound Q < 1.08. The difference An of the spectral indices of the CMB and 
the matter power spectra gives rise to another constraint on the model parameter p and the coupling Q. 

Figure [5] illustrates the bounds derived from the conditions An < 0.05 and s < 2, as well as those from local gravity 
constraints. The models with p close to 1 satisfy all these requirements. It will be certainly of interest to place more 
stringent constraints on the values of p and Q by using the recent data of matter power spectrum, CMB and lyman 
alpha forest. Moreover, the future survey of weak lensing may find some evidence of an anisotropic stress between 
gravitational potentials $ and 'f, which can be a powerful tool to distinguish modified gravity models from the ACDM 
cosmology. 
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APPENDIX A: STABILITY ANALYSIS 



In this Appendix, we briefly summarizes the results of the stability analysis which are necessary to obtain the 
background cosmological dynamics discussed in Sec. IIIII 



1. Stability of the fixed points 

When A is a constant, one can analyze the stability of the critical points (a)-(e) [i.e., Eqs. (f3"2"| - (|37[) ] by considering 
small perturbations 5x\ and Sx 2 around them (2j. We write the equations for perturbations in the form 



— ( s x Xi )=m( f 1 

dN \Sx 2 J \ 5x 2 

and derive eigenvalues /ii and fi 2 of the matrix M. to assess the stability of fixed points. They are given by 
• (a) 



(bl) 



(b2) 



(c) 



(d) 



(e) 



Mi 



3-2Q 2 
2(1 -2Q 2 ) 



M2 



3 + 2QA-6Q 2 
2(1 -2Q 2 ) 



Mi 



3(V6 + 4Q-A) 3 + VQQ 



Mi 



V6 + 6Q 

3(V6-4Q + A) 
V6-6Q 



/'2 



M2 = 



1 + V6Q 

3- V6Q 
1-V6Q 



Mi = 



6 - A 2 + 8QA - ieg 2 
2(1 - 4Q 2 + QX) 



M2 



3 - A 2 + 7QX - 12Q 2 
1 - 4Q 2 + QX 



Ml, 2 



3(2Q - A) 
4A 



1±4/1 



2 - 2QX - 3)(12Q 2 + A 2 - 7QA - 3) 
3(2Q - A) 2 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



Mi = M2 = -3. 



(A7) 
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2. Stability of the de-Sitter point for the variable A 

While the point (e) is stable for constant A, it is not obvious that this property also holds for a varying A. In what 
follows we shall discuss the stability of the de-Sitter point. 

It is convenient to consider the variable \(<j>) as a function of F((f>), i.e., A = X(F). We define a variable, X4 = F, 
that satisfies the following equation 

^ = -2V6Qx lX4 , (A8) 

where the r.h.s. vanishes at the de-Sitter point (e). Considering the 3x3 matrix for perturbations Sx\, 8x2 and Sx^ 
around the point (e), we obtain the eigenvalues 



Mi = -3, M2.3 = -2 



(A9) 



where, Ft = F(4>i) is the value of F at the de-Sitter point with the field value <j>i. Since, Ft > 0, we find that the 
de-Sitter point is stable for 

Q^(fi)>0, i.e., ^(<M<0. (A10) 

We checked that this agrees with the stability condition derived in Refs. [51| by considering metric perturbations 
about the de-Sitter point. 

In f(R) gravity this condition translates into dX/dF < 0. Since in this case, F = e 2 */^ = df/dR and V = 
(RF — /)/2, we have A = —Rf,R/V&V. Then, together with the fact that Rf,R = 2/ holds for the de-Sitter point, 
the condition, dX/dF < 0, is equivalent to R < f,n/f,RR. For positive R this gives 

<^<l, (All) 
J.R 

which agrees with the stability condition for the de-Sitter point derived in Ref. [Tlj . 
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